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1. Introduction 



The measurement of the high mass of M = 1.97 ±0.04 M© for the pulsar PSR J1614-2230 
||l|] has reinforced the question whether the phase transition to hyperon matter |Q or to deconfined 
quark matter [^] in compact stars can be excluded (see, e.g., |5|, 0, ^, IToj, [III, l^^)- 

In the present contribution, we want to discuss an answer to this question on the basis of a 
recently developed approach to quark matter within a nonlocal Nambu-Jona-Lasinio (NJL) model 



coupled to the Polyakov loop [ ]13| , |14| , |15P (denoted here as "nl-PNJL") in its extension including 
wave function renormalization (WFR) [|l^, |l7, 18] which goes far beyond the local PNJL models. 

These models have been developed and tested against modem LQCD data in the finite tem- 
perature domain and at small chemical potentials where also results from LQCD within Taylor 
expansion techniques exist [|l^]. The extension of nl-PNJL models to the full QCD phase diagram 
poses a challenge pH, El]]. In particular their extension to the zero temperature case for studies of 



deconfinement in compact stars has so far been performed without WFR [22, 23]. The inclusion of 
the latter into these studies is a novelty on which we present first results in this contribution. 



2. Nonlocal PNJL model 



The Lagrangian for the SU (2) j nonlocal models, including vector channel interactions, is 
given by 

^ = q{iII)-mo)q + ^M + '^{^) . (2.1) 

where q is the Nf = 2 fermion doublet q = {u,dY , and wiq is the current quark mass (we consider 
isospin symmetry, that is mo = m„ = md). The covariant derivative is defined as = d^— iA^, 
where are color gauge fields. The nonlocal interaction channels are given by 



-Z-mt — T- 



js{x)js{x) + jp{x)jp{x) - ip{x)ip{ 



-^jv{x)jv{x) 



where the nonlocal currents are 



ja{x) = J d^Z g{z) ^ (■^+ I) r„ ^ 

jp{x) = I d'^z f{z) ^ (■^+ I) ^ 9 



z 

X 

2 



a = S,P,V 



(2.2) 



(2.3) 



with Fa = (r5,rp,ry) = (l,/75f , 7o) for scalar, pseudoscalar and vector currents respectively, and 
u{x')dv{x) 



■ u{x')dxv{x) — dxiu{x')v{x). The functions g{z) and f{z) in (23) are nonlocal covariant 
form factors characterizing the corresponding interactions. The scalar-isoscalar current js{x) will 
generate the momentum dependent quark mass in the quark propagator, while the "momentum" 
current, jp{x), will be responsible for a momentum dependent WFR of the propagator. Note that the 
relative strength between both interaction temis is controlled by the mass parameter Kp introduced 



in (2.3) 



In what follows it is convenient to Fourier transform into momentum space. Since we are 
interested in studying the chiral phase transition we extend the effective action to finite temperature 
T and quark chemical potential jj. using the Matsubara formalism. 
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Figure 1: (Color online) Left: Normalized dynamical masses and wave function renormalization for the 
choice of form factors used in this study (set C), fitted to lattice data of Ref. |28| according to p?]]. For 
comparison, more recent data [ p9| ] and alternative formfactors according to set B of Ref. [|l6| are shown. 
Right: The phase diagram of the nl-PNJL in meanfield approximation for the formfactor parametrization of 
set C, see also Ref. |]20|]. Reproducing the slope of the pseudocritical temperature fc = 0.059 from lattice 
QCD for small chemical potentials [ [l9| ] requires a scaled vector coupling Tj — 0.555. 



Concerning the gluon fields we use tlie same prescription as in previous works [14, 16], but in 



our present case we have chosen a /i -dependent logarithmic effective potential described in [ |24| ] 



■30^), 



(2.4) 



where ao = —1.85, a\ = — 1.44 x 10^ , 02 = —0.08, = —0.40. In the present work we set the Tq 
parameter by using the value corresponding to two flavours Tq = 208 MeV, as it has been suggested 
in [H], and used in subsequent approaches, including the nonlocal PNJL [26] and Polyakov 



loop-DSE models 

Finally, to fully specify the nonlocal models under consideration we fix the model parameters 
as well as the form factors g{q) and f{q) following [16, |^], from where we choose the more 
realistic combination of Lorentzian functions (set C, according to the notation in Ref. [|^]). These 
form factors have been chosen [^] such as to reproduce the dynamical mass function M{p) and 



the WFR Z{p) from lattice QCD simulations of the quark propagator in the vacuum [g8|], shown in 
the left panel of Fig. [l[ For comparison, also the more recent data [29] and alternative formfactors 



according to set B of Ref. [16] are shown. 



Within this framework the mean field thermodynamical potential Q^^^ results 



MFA 



-4r£ 



In 



2Gs 



(0^ 

2G^ 



{2ny 

where M{p) and Z{p) are given by 

M{p) = Z{p) [m + ai g{p)] , Z{p) = [1 - 02 f{p)r' 



+ ^{<^,T), (2.5) 



(2.6) 
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In addition, as in [16], we have considered 

, 2 



{2n+\)KT -i}X + <^c 



\p\ (2.7) 



where the quantities (pc are given by the relation = diag(0r, 0^, 0^). Namely, (j)^ = c 03 with 
c = 1, —1,0 for r,g, ft, respectively. In the case of (^p^pj we have used the same definition as in 
but shifting the chemical potential according to 

fL = H - (0 g{p) Z{p) . (2.8) 



The expression for D?^^^ as given in ( 2.5 ) turns out to be divergent and, thus, needs to be regu- 
larized. For this purpose we use the same prescription as in [jl^, A necessary condition to find 
the set of mean field values for ai,2, (0 and 03 which would correspond to the absolute minimum 
of the regularized thermodynamic potential is the fulfillment of the coupled "gap" equations 

anj^f g^MFA g^MFA g^MFA 

doi da2 dco dips 

Two of these meanfields are regarded as order parameters as their nonvanishing value signals the 
breaking of a symmetry in the system: the quark mass gap Oi stands for the breaking of the chiral 
symmetry and the traced Polyakov loop <I> = [1 + 2cos(03/r)]/3 signals the SU(3) center symme- 
try breaking signalUng deconfinement. In the right panel of Fig. [l] we show the phase diagram of 
the chiral transition in the T — /i plane for two values of the scaled vector coupling rj = Gy/Gs- 
The dotted lines show the pseudocritical temperatures of the chiral crossover transition which at the 
critical endpoints (CEP), shown as filled circles, go over to the critical temperatures (dashed and 
solid lines) of the first order transition where the chiral order parameter is discontinuous. Within 
the set C parametrization the value Tj = 0.555 is required for reproducing the slope of the pseud- 
ocritical temperature K = 0.059 obtained in LQCD simulations for small chemical potentials [|l^]. 
In comparison to a vanishing, or lower, vector coupling the CEP is situated at a lower temperature 
and higher chemical potential. The CEP position also depends on the parametrization of the PNJL 



model. In the local limit it even vanishes, see Ref. [gOp for details. 
The EoS for quark matter is then obtained from 



P{T,n) = - min ^ D^^\T,n;{auO2,(0,<t>3}) , (2.10) 



which holds for a homogeneous system. On the left panel of Fig. g we show the 7 = EoS of the 
nl-PNJL model in set C parametrization for different values of tj . 



3. Deconfinement phase transition in compact stars 

The deconfinement phase transition between nuclear matter and quark matter under compact 
star conditions' is constructed by solving the corresponding condition for phase equilibrium 

Pgi^c) = Ph{^c) , (3.1) 

'Here we consider local charge neutrality and /3— equilibrium. A more elaborated treatment considers finite size 
structures ("pasta") in the mixture of hadronic and quark matter. See pip for a recent discussion and references. 
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Figure 2: (Color online) Left: Pressure vs. baryochemical potential EoS for neutron star matter in the 
hadronic phase (DBHF with and without excluded volume Vex, APR) compared to that of the quark matter 
phase (nl-PNJL) for different values of Tj. For Tj ^ 0.5 nuclear matter is stable at low densities. For a phase 
transition to quark matter at high densities (chemical potentials) the vector coupling has to get reduced to 
a value ?]> ^ 0.3. Right: Hybrid EoS, where the quark sector is provided using the interpolation (3.2) and 
two alternatives for hadronic phase are used (APR and DBHF). The phase transition from hadronic to quark 
matter is determined by Eq. (3.1). The green hatched area corresponds to the EoS constraint for symmetric 
matter from the flow data analysis [^sl] shifted by the universal symmetry energy contribution 



where /if is the critical chemical potential for the transition. The index H stands for "hadronic" and 



here we will use one of the standard nuclear matter EoS: APR [g2|] or DBHF \^%. The index Q 
stands for "quark" and the corresponding EoS will be introduced based on the following discussion. 
Inspecting the left panel of Fig. |^ we observe that choosing for Pq{^) = P{0,IJ.), i.e. the EoS 
( |2.10 ) for the nl-PNJL with a constant vector coupling strength rj there occurs one of the following 
problems: (a) quark matter is favorable for all chemical potentials, even at lowest densities (for 
rj ~ 0.3); (b) nuclear matter is favorable for all chemical potentials, even at highest densities (for 
77 ~ 0.5); (c) quark and hadronic EoS cross each other either in the "wrong" way (confinement at 
high densities) or more than once (for rj ~ 0.4 . . . 0.5) since they are marginally indistinguishable 
in the relevant range of chemical potentials. The latter situation has been called the "masquerade" 



problem []34[]. 

In this situation we suggest the following procedure to avoid unphysical solutions and the 
masquerade problem. We anticipate a chemical potential dependence of the scaled vector coupling, 
77 (/i), in such a way that at low chemical potentials the value 77< = 0.555 is set which is necessary 
to reproduce the slope of the pseudocritical temperature and at high chemical potentials, beyond 
a fiducial range for the deconfinement transition, an asymptotic value is attained which makes the 
quark matter EoS favorable over the hadronic one under consideration. For H = DBHF, this is 
77> < 0.3. For the transition between both regions we suggest here an interpolation ansatz 



Pe(At) = P(0,At;77<)/<(At)+mM;'7>)/>(M) 



/<(M) 



1 =Ftanh 



(3.2) 
(3.3) 



For the two free parameters of the interpolation we suggest the ranges /i '-^ 360. . .410 MeV and 
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Figure 3: (Color online) Hybrid star sequences based on the DBHF and APR EoS for nuclear matter and 
the nl-PNJL model for quark matter proposed here, shown together with constraints on the maximum mass 
from PSR J1614-2230 [III] and on the mass-radius relation from RX J1856.5-3754 [^. 



r ~ 30 MeV. The motivation of this choice comes from considering the influence of a nucleonic 
excluded volume, shown also in the left panel of Fig. which stems from the quark Pauli blocking 
between nucleons and thus estimates the onset of quark substructure effects like the overlap of 
hadronic wave functions in dense matter modifying effective quark-meson couplings. 

Two sets of resulting hybrid EoS, based on DBHF and APR as hadronic EoS, respectively, are 
obtained by the Maxwell construction according to Eq. (3.1). They are shown on the right panel 
of Fig. ^ together with a constraint for compact star EoS obtained from the flow data analysis for 
symmetric matter [35| shifted by a contribution due to the universal symmetry energy of Ref. [^. 

With these hybrid EoS the Tolman-Oppenheimer-Volkoff equations [37, 38] have been solved 
and the resulting sequences of nonrotating compact star configurations have been obtained. The 
corresponding dependencies of the mass M of the stars on their radius R and on their central density 
n{0) are shown the left and right panels of Fig. ^, respectively. 



The hybrid star sequences based on the APR EoS do not fulfill the constraint on the maximum 
mass from the mass measured for PSR J1614-2230 and also do not obey the constraint on the 



minimal radius for a given mass from RX J1856. 5-3754 [ |39| ] since these hybrid EoS describe too 
"soft" compact star matter ^. The hybrid stars based on the DBHF EoS fulfill both constraints 
and, depending on the specific /i— dependence of the vector coupling strength T](/x), can have an 
onset of quark matter in their interiors already at a mass as low as 1.2 Mq. This would include 



^Note that the analysis of X-ray burst sources in Ref. can provide at best a lower limit on the star radius which 
is embedded in that for RX J1856. 5-3754, see Ref. [ pi] ] for a discussion. 
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all compact stars with measured masses known so far. However, these EoS are too stiff for the 



constraint from flow data [ ]35| ] combined with the universal symmetry energy contribution [|36|]. 



4. Results and Conclusions 

A nonlocal Polyakov-NJL model with vector meson meanfield has been calibrated with lattice 
QCD data for the momentum dependence of the wave function renormalization and dynamical 
mass of the quark propagator in vacuum as well as for the slope of the pseudocritical temperature 
at small chemical potentials which fixes the vector coupling strength under these conditions. 

A straightforward application to nonzero temperatures and chemical potentials within the Mat- 
subara formalism results in a prediction for the position of the critical endpoint. 

Considering compact star matter at T = results in the conclusion that quark matter in com- 
pact stars cannot occur since its pressure stays always below that of standard hadronic EoS for 
all relevant baryon densities, unless a medium dependence of basic nl-PNJL model parameters is 
invoked. A generic ansatz for the chemical potential dependence of the vector coupling strength is 
presented which allows to estimate the ballpark for such a dependence allowing hybrid star con- 
figurations which fulfill basic compact star constraints. Examples illustrate the sensitivity to the 
nuclear EoS: While the set of hybrid star EoS based on APR is too soft to fulfill mass and radius 
constraints, those based on DBHF satisfy these constraints but appear too stiff for the flow data 
from heavy-ion collisions. Thus compact star observations and heavy-ion collisions together may 
guide the path towards the EoS of dense QCD matter. 
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